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INTRODUCTION 

In .1962 Julian Schwingor showed Ij that the photon of the 1+1 dimensional Quantum Elec- 
Ramies (QED) ac q o,es * one loop a m ass ™, -/Vi- —nee of the si„ gul a r 
infrared behavior of the fermionic propagator. A little later, Anderson |2j showed that photons in 
conductive media may propagate as massive excitations. On the other hand, the gap generation 
in superconductors has as a consequence the Meissner effect, according to which static magnetic 
fields are exponentially suppressed within superconductors. Phenomenologically, the gap can be 
described by the scalar field condensate of the Landau-Ginzburg model, where the role of the scalar 
is taken by Cooper pairs. Abelian and non-Abelian gauge theories in thermal equilibrium have 
also been studied extensively It has been found that, while electric fields are Debye screened, 
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and cannot freely propagate, magnetic fields (at the perturbative level) do propagate as massless 
excitations, and screening is only dynamical. A realization that scalar condensates may be built 
by a fundamental scalar field, led to the establishment of the celebrated Higgs mechanism 0, j], |j| 
(sometimes referred to as the Anderson-Englert-Brout-Higgs-Kibble mechanism), which currently 
represents the most popular explanation for the origin of the mass for the electroweak gauge bosons 
and fermions. This is of course not the whole story. A large part of the mass of mesons and baryons 
can be attributed to nonperturbative effects of a strongly coupled QCD. For an insightful histor- 
ical overview of the development of our understanding of the concept of mass we refer to |7| ■ We 
paraphrase Okun: A tiny photon mass, albeit gauge non-invariant, does not destroy the renormal- 
izability of QED j], Q], and its presence would not spoil the beautiful agreement between QED 
and experiment. This motivates (so far unsuccessful) incessant searches for a nonvanishing photon 
mass |ld, llll]. 
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aite recently, a novel mechanism for mass generation of gauge fields has been discovered 
, which is operative in the presence of rapidly (superluminally) expanding background space- 
times. The concrete model within which the mechanism has been studied is scalar electrodynamics 
(<5QED), endowed with a massless, minimally coupled, scalar field evolving in a locally de Sitter 
space-time background. The photon mass is radiatively induced (at one loop) as a consequence 
of the canonical photon coupling to the amplified charged scalar fluctuations, which is in contrast 
to the Higgs mechanism, where the mass is induced by scalar condensate. We now compare this 
gravity induced mechanism with the case of thermal media. While magnetic field s ge t initially 
suppressed (with respect to the conformal vacuum), at asymptotically late times Jjj|, just like 
in thermal equilibrium, they are not screened, and the final amplitude remains comparable to 
the conformal vacuum, B ~ B va , c oc a~ 2 , where a denotes the scale factor. On the other hand, 
the photon mass in a rapidly expanding Universe implies exponentially enhanced (anti- screened) 
electric fields, E oc a~ 3 / 2 , which is in contrast to the Debye-screened electric fields of thermal 
media. 

,0 



Here we reconsider the mass generation mechanism of Refs. 
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14, UM in the context of 



massive scalar electrodynamics, whose Lagrangean density in a general metric field reads, 

£$qed = - l -^gg» p g u °F^F pu - V^gg^(D^D u <p - ^g(m\ + £i?)<f^ , (1) 



where = d fJi A u — d u A fJi is the field strength, D« = + ieA^ the covariant derivative, 
det[<T,], a nd R is the curvature scalar. Of course, through the Coleman- Weinberg mechanism 
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101. 



the effective <1>QED acquires radiative corrections, which can be expressed through a renormalized 
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mass — > tur and a renormalized quartic coupling, A#. Since we shall not study the consequences 
of radiatively induced symmetry breaking in the scalar sector, for our purposes it suffices to consider 
the physical effects of a nonvanishing scalar mass, > 0. We do however discuss briefly what 
happens in the case when m 2 . < (or, more precisely, when + £i? < 0). 
When recast in the locally de Sitter background, 

= a 2 (r))r]^ , rj^ = diag(-l, 1, 1, 1) , (2) 

where a denotes the scale factor, which - when written in terms of conformal time r] - reads 

«W = "^ fa<0), (3) 

£$>qed in Eq. (^Q) reduces to, 

Uqed — > —aPS^F^Fpa ~ a D - 2 ^ v {d^ - ieAJtf {d v + ieA v )<t> - a D (m 2 + £R)<p*0 , (4) 

where H is the Hubble parameter in inflation, t?„ = (cL,V), D denotes the dimensionality of 
space-time, and we made use of y/—g — ► a D . Recall that in D-dimensional de Sitter space-time, 
R = D(D — 1)H 2 . In this paper we shall assume a nearly minimally coupled, light scalar, for 
which ^ < 1 and <C H 2 . Prom Eq. © it immediately follows that, in D = 4, gauge fields are 
conformally invariant, while conformal invariance of scalar fields is granted for £ = 1/6. 

In section[n]we introduce the scalar two-point function G = G(y), which (in D-dimensions and 
in the absence of electromagnetic fields) can be expressed in terms of the de Sitter invariant length, 
y, as the following hyper geometric function, 

<;•-' . „) r (^l_„) „(D-1 D-l D y 



where v is given in Eq. (|2U|) . y = y(x;x') can be expressed in terms of the conformal coordinate 
interval, Ax 2 (x;x') = rf l '{x il — x')(x v — x' u ) = —{rj — i]') 2 + ||x — x'\\ 2 , and the geodesic distance 
£ = £(x;x'), 

y(x; x') = aa'H 2 Ax 2 (x;x')=4 sin 2 (^H£(x; x') V (6) 

where a = o(ry), a' = a(rf). 

It is instructive to expand the solution © around the minimally coupled case, such that in 
D = 4 we get, 
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where the parameter |s| <S 1 is defined by, 



D - 1 



(g-i )g3 +Q (W + ^)/ g 



(8) 



H 2 

This is to be compared with the corresponding expression for the massless two-point function in 
D = 4, 

Apart from a constant term, the two solutions (j7J) and © differ by the term \H 2 / (87 T 2 )J ln( aa' ) 
which is responsible for breakdown of the de Sitter invariance by the massless propagator Il2.ll5l. ll7l . 
llsUlfit l2(il2l|. and it is a consequence of the growth of scalar fluctuations during inflation 22, 23]. 
Note that, in the limit when s — > 0, the propagator (pj) becomes formally singular, and hence 
ill defined. This is a consequence of the fact that it is not possible to construct a nontrivial de 
Sitter invariant propagator for a massless scalar 17J, |2J]|. Note also that, since the non-Hadamard 



terms in the propagators ® and © dominate in the infrared, they play a similar role in our 
mass generation mechanism as the fermionic states in the Schwinger mechanism, whose infrared 
behavior is sufficiently singular to lead to mass generation in 1 + 1 dimensions. 

In order to investigate how is the mass generation mechanism affected by the difference in the 
scalar propagators, in section ITTTI we outline our calculation of the one-loop, renormalized, vacuum 
polarization tensor of <£QED to order s°. We then use standard techniques to arrive at the retarded 
tensor, which is then used in section Hvl to study the photon dynamics. 

Our main result is presented in section HVl and can be summarized as follows. To order s°, the 
photon dynamics during inflation is governed by the Proca Lagrangean, which in D=4 reads, 

£p r oca = —r^rTF^F^ - ^a 2 m 2 y ' A^A V + 0(s°) (10) 

with the photon mass-squared (accurate to order s°), 

2 aH 2 _ 3aH 4 
m ^~^T~ vr(m2+6R) ' ( } 

where a = e 2 /4ir is the fine structure constant. The Proca Lagrangean is derived in a generalized 
Lorentz gauge, 

3 M (aVM*) =0. (12) 

This completely fixes the gauge, such that the photon contains two transverse and one longitudinal 
degree of freedom. In analogy to the Higgs mechanism, the longitudinal degree is supplied by the 
charged scalar fluctuations. 
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A consequence of the singular behavior of the two-point function (JJJ) in the limit when s — > 

12 , Q] is then to use 



is an ill-defined photon mass The correct treatment of the problem 



the massless propagator @, in which case one recovers the (finite) result of Ref. [la], which can 
be interpreted as the (electric) photon mass, m 2 = (2a / ir) H 2 ln(a) . The time dependence is a 
consequence of the broken de Sitter invariance by the massless scalar propagator Q, which is 0(4) 
invariant. We expect that we would get a similar behaviour for a light nearly minimally coupled 
scalar, were we using a scalar propagator that is invariant only under the subgroup 0(4) of the de 
Sitter group SO(A, 1). The divergent behaviour of the photon mass found in Ref. [la] in the limit 
when a — > co is related to the m \ + £JZ — > behaviour of the photon mass 1)1 ljl in the de Sitter 
invariant case. 

Furthermore, it is interesting to compare the cases when s is small but positive, with the case 
when s is small but negative. In the case 1 3> s > 0, the photon acquires a mass and 
the photon amplitude oscillates and gets suppressed during inflation. On the other hand, when 
— 1 <C s < 0, the mass-squared 1)11(1 becomes negative, indicating instability associated with a 
growth of the photon wave function. Consequently, when compared with the conformal vacuum, 
both electric and magnetic field can grow during inflation. This instability may imply primordial 
magnetic field generation of large amplitudes and on cosmological scales. The quantitative details of 
this mechanism will be addressed elsewhere. Another crucial difference between the mass generation 
in the case of a minimally coupled scalar and that of a near minimally coupled, light scalar is in 
the following. While in both cases the electric field is enhanced (anti-screened) with respect to 
the conformal vacuum, the magnetic field dynamics is strongly affected only in the latter case, in 
which the field is screened as B oc a^^+jV^K/^) 2 _ 

Another noteworthy consequence of the photon mass generation is implied by considering the 
photon dispersion relation, which can be extracted from the Proca equation (obtained by the 
variation of Eq. ()10() ). 

rfPd p {d v A^ - d^A u ) - a 2 m 2 A^ = . (13) 

Assuming a dependence e tk ' x on the spatial coordinates and decomposing the photon field into 
transverse and longitudinal components, one obtains different equations of motion for these modes. 
However, in adiabatic limit, discussed in section IIV A( which is valid for non-relativistic modes 
(fcph = [|fc[|/o <C m 7 ) if m 7 3> H and for relativistic modes (A: p h 3> m 7 ) if k p ^ 3> H (for longitudinal 
modes) or 3> (Hm^) 1 / 2 (for transverse modes), one obtains the same dispersion relation 



u ~ y k 2 + a 2 m 2 (14) 
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for transverse and longitudinal modes. The propagation speed of massive photons is then given by 
the group velocity 



(15) 



which, for ||A: p h|| <C m~ can be < 1 (recall that in our units, the speed of light in vacuo, c = 1). 
Here we used the standard definitions for the physical frequency, u; p h — u/a, and the physical wave 
vector, kph = k/a. Since, in the case when < s « 1, m 7 can be S> H, a large class of both 
sub- and superhorizon photon modes may propagate with speeds, which are much smaller than the 
speed of light in vacuo. We have thus discovered that, provided inflation lasts a sufficiently long 
time, such that the de Sitter invariant solution for the scalar two-point function has had enough 
time to get established, light can propagate very slowly during inflation. Moreover, because the 
physical momentum scales as k p h oc k/a = ke~ N , where N = Ht denotes the number of e- foldings, 
t the physical (cosmological) time, and k the conformal momentum (the physical momentum at 
time t = 0, or rj = —1/H), the group velocity at asymptotically late times drops exponentially with 
time, Wgroup oc e~ Ht (t — ► oo). We also note that, if the adiabatic conditions are not met, one can 
still calculate u gI . oup = dio/dk by solving the equation for u = u;(k,rj) exactly. As an example, in 
section ITV Al we show that the relativistic longitudinal photons, for which k p h,H 3> m 7 , propagate 
sup er luminally. 

The mass-induced slow-down of light differs significantly from the light slow-down reported not 
a long time ago 27] in ultra low temperature sodium gas, where the slow-down is induced by a 
steep frequency dependence of the index of refraction. Since the effect is of a resonant origin, 
it pertains only in a very narrow range of frequencies, which is in contrast to the mass-induced 
slow-down during inflation, which is effective for a broad range of frequencies. We also mention the 
work of Ref. 30], where an extreme version of the Fresnel effect is considered, in which superfluid 
vortices moving at a superluminal speed are used to generate an optical Aharonov-Bohm effect. 
Moreover, if they rotate faster than the speed of light in the medium, they can trap the slow light, 
and in this sense mimic the black holes of general relativity. 

Finally, in section ^ we present an estimate of the cosmological magnetic fields produced as a 
consequence of a radiatively induced photon mass during inflation. Our conclusion is that the field 
strengths thus produced are sufficiently strong to satisfy the bounds on the seed magnetic field of 
the galactic dynamo mechanism. Appendices are reserved for technical details of the calculations. 
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II. SCALAR PROPAGATOR 



A complex scalar field of scalar electrodynamics (@J) obeys (in the absence of electromagnetic 
fields) the following equation of motion, 



rTd^a D - 2 d, 



a D (m 2 



£R)(f) = 0. 



(16) 



Of course, the two-point Wightman functions, (a\4>(x)4>^ (x')\a) and (a\<t>' (x')(j)(x)\a) , satisfy the 
same differential equation. For a de Sitter invariant vacuum state la), where a represents the (rea l) 



parameter that classifies all globally de Sitter invariant vacua |2J| (c/. also Refs. |3l| and 32]), 
they can depend only on the geodesic distance between x and x', and may be written as a function 
of z(x,x') = 1 — y(x,x')/4 only (c/. Eq. ©)• For such a function G(z) the differential equation 
can be recast as, 



rTid^jduz) d 2 
a 2 dz 2 



G + 



rf v d^ v z (D - 2)H(d z) 



d 

dz 



G - (m 2 + £R)G = . 



(17) 



cr a 

Upon calculating the derivatives of z, and contracting with rj^ 11 , one obtains the hypergeometric 
differential equation, 

d 



2(1 _ 2) ^G + D (l_ 2j _ G 



(mi+{fl)G = 0. 



(18) 



G(y) = G(z) in terms of the hypergeometric functions 

D 



2J.I25, 26], 



G{y) =c 2 Fx( + v, ■ 



u, — ; 1 + c o^i h v, 

2 2 2 4/ V 2 2 

where (so far) c and d are arbitrary constants (which are dependent on a), and 



,.f;f), m 



D - 1 



ml + gi? 
H 2 



(20) 



The hypergeometric function becomes singular |2l(, when the last argument approaches 1 or 
— 1, and has a branch cut from 1 (—1) to +oo (— oo). In order to completely specify the Wightman 
functions, we demand that there be a singularity only if x and x' are light-like related. Furthermore, 
at short distances this singularity should have the Hadamard form of the Minkowski space two-point 
functions, since (at short distances) the scalar field can be only weakly affected by the expansion 



of space-time. With these assumptions one finds d = 0, such that ()19|) reduces to 

1 



G(y) 



(47r)Tr(f) 



^h d - 2 2Fl 



D-l D 



D 
v, — ; 1 

' 2 ' 



(21) 
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where the precise value of c is dictated by the canonical commutation relation of (f> and its canonical 
momentum j^jj. This choice corresponds to the unique vacuum, \a = — oo). (The a = — oo vacuum 
is in literature known under various names: Chernikov-Tagirov vacuum, Bunch-Davies vacuum, 
Euclidean vacuum, thermal vacuum.) The vacua with a ^ — oo contain an additional singularity 
at the antipodal point, y = 4 (which corresponds to rj = —77 and x' = x), which can be also 
interpreted as a constant (a-dependent) nonvanishing particle number at all momenta, such that 
a / -00 vacua are associated to a divergent stress-energy tensor. While this divergence should be 
regularized, it can be done so only for a particular choice of a. In our point of view, the requirement 
that, at short distances, one should recover the Hadamard form for the singularity of two-point 
functions (which, by the way, has been tested in particle accelerators), singles out the a = —00 
vacuum in a natural way. Indeed, the Hadamard form can be violated only at the price of giving 
up locality in position space. The nonlocality is namely necessary if the Hubble scale physics is to 
influence the physics at much shorter scales, responsible for the Hadamard singularity. 

Analogous to different pole prescriptions in the momentum space integrals for the various forms 



of the scalar propagator 



441 ] . one can write the (anti-)Feynman and Wightman propagators as the 



same function G of the appropriately modified de Sitter length functions, 



iA w (x,x') = G(y w ) (b,b' = +,-), 



(22) 



where 



iA(x,x') = iA ++ (x,x') = (0\T[cj)(x)$ (x' 
iA + _(x,x') = (0^ (x')(f)(x)\0) 
iA_ + (x,x') = (O|0(a?)0 t (aj')|O) 
iA__(a;,x / ) = (0|f [0(x)0 f (x' 

Here T (T) denotes time (anti-time) ordering and 



yw 



1 i 

TjTj 



with 



(Feynman) 
(Wightman) 
(Wightman) 
(ant i— Feynman) 



Ax 2 = Ax 2 



Axi 



Ax' 



Ax 1 



(| 77 — t/| — id) 2 + ||x — cc'|| 2 , 
(77 — 77' + i5) 2 + \\x — x\\ 2 , 
(77 — 77' — i5) 2 + \\x — x\\ 2 , 
(\v ~ v'\ + + \\% ~~ ^'l| 2 - 



(23) 
(24) 
(25) 
(26) 

(27) 

(28) 
(29) 
(30) 
(31) 
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From now on we shall use the notation, 

Ax 2 , s 

y = y++ = — • 32 

r/r/' 

In order to make progress toward calculating the vacuum polarization tensor in dimensional 
regularization, we shall expand the scalar propagator (|21 j) — (|26|) around the massless minimally 
coupled case and in the vicinity of D = 4. To this purpose, it is useful to define the parameter s 
as, 

s = — v, (33) 

where v is defined in (|20j) . and the parameter 

e = 4 - D . (34) 

For a light and/or nearly minimally coupled scalar, for which <C H and |£| <C 1, we can expand 
v in powers of (mi + £,R)/H 2 , to obtain the following approximation for s, 

In AppendixEJwe show (see Eq. (|A8|) ) that, when expanded in s, y and e, the Feynman propaga- 
tor (|22H23|) reduces to, 

iA(x, x') = 0f(y) + s-^^ [g(y) + %)] + 0(e, s 2 ) , (36) 



where 

P 



e\yi 2 e r(3-e 



1 /7T£ 



r(2-f 

2 - n(n + 2) ln(y/4) /y\ n 



n=l 

Here T denotes the Euler-Gamma function, 7e ~ 0.577 the Euler-Mascheroni constant, ip is defined 
by tp(z) = V (z) /T(z) and C(s,e) is an order s term that is independent on the coordinates. We 
will later find that it does not contribute to the vacuum polarization at order s°. The 0(e) term in 
the infinite sum is not important, because it will not be needed for the dimensional regularization 
and renormalization of the vacuum polarization tensor outlined in section ITTT1 and Appendix IO 
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For calculating the one-loop vacuum polarization tensor in a locally de Sitter space-time, we 
will also need the coincident limit x' — * x (equivalently y — * 0) of the propagator (|3fij) . To obtain 
this limit in dimensional regularization, one must assume that e is large enough that a vanishing 
quantity like y is raised to non-negative powers only. Thus we get 

Jim iA(x, x') = P ^Z&) [\ + * cot (?) " 7E " ^ (3 " £) ] + ° (S) • (41) 

This is in contrast with the massless propagator @, which grows logarithmically with the scale 
factor during inflation, iA(x,x) m _>o °c m(a). 

III. ONE-LOOP VACUUM POLARIZATION IN A LOCALLY DE SITTER SPACE 

We are now ready to calculate the vacuum polarization tensor of scalar electrodynamics (^QED) 
in curved space-time backgrounds (^Q), @ in the one-loo p ap proximation, based on which one can 
study the dynamics of photons during inflation (c/. Ref. [l2j]). The relevant diagrams contributing 
to the one-loop polarization tensor are shown in figure ^ The counterterm in figure ^ (3) is added 
in order to cancel the divergence, appearing in dimensional regularization in the limit when D — > 4. 
Using the position-space Feynman rules of Appendix El one finds the following expressions for the 
one-loop graphs shown in figure Q 

i[^U u ]^(x,x') = -2ie 2 ^/^gJx)g^{x)iA{x,x)5 D {x-x'), (42) 
= - 2ie 2 a D ~ A(s, x)8 D (x - x') (43) 

i[^f\x,x') = 2e 2 y^gJx)g^(x)y^gJx T )g^(x') 

x [(d p iA(x, x'))d' a iA{x, x') - iA{x, x')d p d'jA(x, x')] (44) 
= 2e 2 a D - 2 a' D - 2 ^ p if a [{d p iA(x, x'))(d'jA(x, x')) - iA(x, x')d p d'jA(x, x')] , (45) 

where iA(x, x') = iA(x', x) is used and 

^lr] (3) (*,*') = - i5Zd p {^gJx)[g^{x)g^{x) - g^(x)g l 'P(x)}d'j D (x - x')), (46) 
= -i5Z[»P v ]a D - i 5 D {x-x'). (47) 

Here a and a 1 are defined by a = a(ry), a' = a(r)'), and [^P"] is the transverse projector defined in 
Eq. dHJ) below. One can use Eq. dHJ) for iA(x,x) in (H21) and (j3H for iA(x,x') in $o§. Then the 
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0) 



fc^X/V (2) 

/ x' 



(3) 



X X 



FIG. 1: The one-loop Feynman graphs contributing to the vacuum polarization tensor in scalar QED. 
sum of the graphs Fig. ^(1) and Fig. [2(2), ^[ M n^ +2 ] can be derived, one obtains (see Appendix IU|) 



*rrrr +2 ] (*,*') 



where 



g r ^(i-f) i 

2vr 3 \ 2(3 -e) L J Ax 4 " 2e L 



2 s 



+ 2 



r/r/'Ax 2 



+ 



HI 



2(y - 4:)r] 2 r]' 2 



+ ^P 1 



4s 2(y-A)J 4 V 4 



+ 0(b) , 



(48) 



[mp-] = v w v ™( Vp<7 d , -d-d' p d a ) 

[tp*] = n^irtijV-V-dldj) 



(49) 



are the manifestly transverse projector operators, ^[^P"] = = ^[^P"], d /i [ /i P 1 '] = = ^[^P"], 
and d' ■ d = rj^d'd^, V • V = c^'Sj (greek indices a, /3, .. run from to 3, while latin indices i, j, .. 
run from 1 to 3). In Eq. (|48|) D was taken to 4 in all terms, except in the first, because only this 
term requires regularization. Ax 2 and y are given by (|28|) and ()32|) . and Lz2 is the dilogarithm 
function defined by L%2{z) = — f? ^ dt. 
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The first term in Eq. (|48|) is exactly what one would find for Fig^(l) and Fig. ^(2) for a 
massless scalar field in flat space [12]. Its divergence for e — > can be seen from 



1 



2 



1 



Ax 4 " 2£ 2e(l-e) Ax 2 ~ 2e 



Combining this with 



4iir 2 2 



Ax 2 - £ T(l - f ) 



6 D {x-x'), 



one finds 



1 



S 1 



Ax 4 ~ 2e 2e(l - e) 



1 



27r 2 i(0F/i)- £ xD 



e(l -e)r(l 



^(x-x'), 



(50) 



(51) 



(52) 



,Ax 2 ~ 2e Ax 2 ~ e 

where [i is a parameter that will be used for regularization. Taking e — > in the first term gives 

2vr 2 f(y?/i)- £ xD 



(P rln( M 2 Ax 2 ^ 



Ax 2 



-S D {x-x') 



(53) 



e(l-e)r(l-f) 

This expression can be used for the first term of (|48|). Its divergent local part can be canceled by 
the counterterm FigQ(3), which can be rewritten as 



^lrpOc, x') = -i8Z[»P v ]{l - eln(o) + 0(e 2 ))5 D (x - x') , 



(54) 



where bZ is chosen such that, to the lowest order in e, the counter-term cancels the divergence in 
Eq. (|53|) . Note that the logarithm in (|54|) still gives a contribution to the renormalized vacuum 
polarization, which grows logarithmically with the scale factor. Thus we find 



a ( 



f ln(fj, 2 Ax 2 ) \ , »tt 4 , 



+ 0(s) 



2s 



+ 2 



r/r/'Ax 2 



+ 



mi; 



2(y - 4)r/ 2 ?/ 2 



4s ' 2(y - 4) y " 4 * 2 V 4 



(55) 



which is the renormalized vacuum polarization tensor for D = A written in a manifestly transverse 
form. 



A. The retarded vacuum polarization tensor 



In this section we use the Schwinger-Keldysh formalism [32 



34, 



35 



to construct the retarded 



vacuum polarization tensor, [ ^IIr4n] , required to study the dynamics of photons during inflation. In 



14 



order to do that, one is led to modify the Feynman rules of Appendix[Bl such that, in addition, the 
vertices become signed as b = + or — , while the propagators (|22H2fi|) (which connect these vertices) 
acquire two signed indices, iAbb'(x,x') (b,b' = +, — ). Repeating the procedure of section ITTT1 in 
which we calculated the Feynman vacuum polarization tensor, i [ ^n^ + ] (x, x') , we shall now outline 
how to derive other relevant vacuum polarization tensors, z[^II^,] . Since the different propagators 
(122|) are the same function of the appropriately modified de Sitter length functions (|27j) . in order 
to get the different versions of the polarization tensor, we just have to use (|2"T|) in A couple 

of subtle remarks are in order. In the derivation of the vacuum polarization tensor we used 

iA(x,x') = iA(x',x). While iA ++ (x, x') and iA (x,x') are symmetric under the exchange of 

x, x' , for the 'off-diagonal' Wightman functions iA^ and iA (_, the following symmetry of the 

propagators ought to be used 



iA w (x, x') = iA Vb {x', x) (b, b' = +, -) 



(56) 



which can be easily established from Ax 2 h ,(x, x') = Axl, b (x', x). Moreover, because the vertices are 
signed, the off-diagonal vacuum polarization tensors, i[^n^_] and i ], acquire an overall minus 

sign with respect to i[^II+ + ]. Finally, there are no H — or — + seagull graphs or countertems, and 
there are no local terms coming from Fig.^(2) [12]. Based on the above considerations and Eq. (|55|) . 
we can write, 

2 \n{ii 2 Ax 2 w ) iir 2 



+ 0(s) 



2 



, + {b + b')Ha)5\x-x') 
24Ax bb , 6 



^ln( yty /4)-l+2 
rirfAx 2 



+ 



ln(y w /4) 
2{y bb > - 4:)rj 2 r]' 2 



l \n 2 ( y -^\+\J y ^' ' 



V 4 



V 4 



1 



+ 



Vw 



4s 2{y w - 4) 



;Li 2 1 



Vw 
4 



(57) 



For the derivation of an effective field equation of the photon field we shall need the retarded 
vacuum polarization tensor, which is defined by 



rn^K*,*') = r"n* + ](x,x') + rn^_](x,x') 



(58) 



Because the two tensors contribute with an overall minus sign (c/. the sign prefactor bb' in Eq. (|57|0 . 
the contributions to the retarded vacuum polarization come only from branch cuts and singularities 
of (fBTI) in Ax 2 h , and yw and from the term oc (b + b') ln(a)5 4 (x — x'). In order to extract these cut 
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and pole contributions, the following formulas are useful, 
-\n(fi 2 Ax 2 bbl ) 



2AAx 2 w 



1 



iln(Y)-^ + 2 
777/ Ax 2 , 



1 



^lnVA^,,) - ^ln(/x 2 Ax^ 



~d\ln\Ax 2 w )) + (3 - ln(4 W ') - i)«9 2 l*(Aa&0 



(59) 



from which it can be seen that the contributions, which are singular at the light-cone, yield finite 
cut and/or pole contributions. The contributions from the logarithms in (|57|) are simply, 



ln(Aa; 2 ) -ln(Ax 2 ) 



In 



2/2/+H 



ln( 



4i7r In 



In 



y^ 



2i7r9(Ar 2 )e(Ar ? ) 



At 2 



ln 2 (Axi + ) -ln 2 (Ax^ 



Uvf e(Ar 2 )e(A,) 
4nrln|Ar 2 | e(Ar 2 )e(A??) , 



(60) 

(61) 
(62) 



where At 2 = Arf — ||Ax|| 2 , O is the Heaviside step function, 0(Ar 2 ) = 0(|Ar/| — ||x||), Arj = Tf — rf 
and Ax = x — x'. From the integral representation of the dilogarithm function one finds, 



Li 2 1 



Moil 



At 



2m In 1 + — — 9(Ar 2 )e(A?7) 



(63) 



4 / V 4 / V 47777' 

Combining these equations we find the following expression for the renormalized, retarded vacuum 
polarization tensor to one-loop order, 

reSJC*. x>) = ^| - [»P«] {±d* [G(Ar 2 )8(Ar ? )(ln |/, 2 Ar 2 | - 1)] + 1 ln(a)5\x - x')} 

e(Ar 2 )e(Ar ? )]} 



i^pui J_( 5 2 \o(A T 2 )e(Ar]) In |A 
4t7t/ I L 

r upv] e(Ar 2 )6(Ar ? ) 

^ (y - 4)r/ 2 ry /2 

. AT 2 ]flf Ar?") r 1 r 

In 



+ [^| 
+ O(s). 



+ (3-ln(4 W / )-^)(^ 



fj*j 2 ^2 



I2 



+ ln(l-| 



+ 



y-2 1 

- + 2 

y-4 2s 



} 



(64) 



IV. EFFECTIVE FIELD EQUATION AND PHOTON MASS 



An effective field equation for photons in de Sitter space-time can be derived by the variation 



of the Schwinger-Keldysh effective action. The result is the following non-local equation 



12], 



7] up r] 



^d u F pu + J d 4 x'[^ n ](x,x')A u (x') + 0(A 3 ) = 0. 



(65) 



For simplicity here we shall approximate the retarded vacuum polarization tensor (|64|) by its leading 
order 0(s _1 ) contribution, 



a 



8vr 2 £ 



„ n 2 e(Ar 2 )9(A?7) 
7/7/ 



„, 29(Ar 2 )e(Ar ? ) 



2 /2 



(66) 
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which we shall use to study the effective field equation (|65|). This is justified provided + £R\ <C 
.ff 2 . Upon neglecting the 0(A 3 ) contributions in (|65|). we seek the solutions in the form, 

^(x') = £,(^r ? , )e^', (67) 

where 

(d' -^)e (k,v') = ik-e(k, V '), (68) 



which is equivalent to the generalized Lorentz gauge (|12[) . This gauge is obtained by requiring that 
divergence of the Proca equation ()69|) in de Sitter space vanishes, as implied by gauge invariance. 
Inserting Eqs. (|66|) and l67|) into (|65|) and evaluating the integral gives the following Proca equation 
(see Appendix El for the derivation) 

^dyidpA^ - d^Ap) - a 2 m 2 A^ = , (69) 

with the photon mass given by 



aH 2 

m 2 = ^r + 0(s 



2 _ 

7TS 

3aH 4 „ff m l + i R \ 



7r(m2+6R) U IP 

where R is the curvature scalar, which in de Sitter space, R = 12H 2 , and £ specifies the coupling 
of the scalar field to gravity. Thus photons that couple to light minimally coupled scalar particles 
acquire a mass in inflation. A remarkable feature of this result is that, even though a <C 1, such 
that the one-loop approximation is justified, the photon mass may be much larger than the Hubble 
parameter. 

A. On the speed of light in inflation 

We can use the Proca equation (|69|) to study the propagation of light in inflation. We seek a 
solution of the form, 

A ll (x)=e ll (k, V )e i ^, (71) 
where e M (fc, rj) = (eo(k,rj),e(k,r])). Thus the fi = component of (|6T))) can be written as 



idnk ■ e(k, ri) 

eo(k,r)) = —g- Yf, (72) 
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which tells us that the zeroth component of the photon field traces the spatial components. Using 
the gauge (|fi8|) and decomposing e into the longitudinal and transverse parts 

e L = ^fi-k , e T = e - e L , (73) 
k 2 

the spatial components of the Proca equation (|69|) can be recast as 

(d 2 + k 2 + m 2 a 2 ) e T = (74) 
d 2 + P + m 2 a 2 + ^ Ha] \ - do) e L = 0. (75) 



k 2 + m%a 



2 p 2 

-r 



Consider now a longitudinally polarised photon Al{x) = £l& x , el = £h{k, rj). While the magnetic 
field is trivially equal to zero, Bl = a _2 V x Al = 0, the electric field does not vanish, 

^ = ~-2 T-JT, V2 • (76) 

a 2 1 + (k/am^) 2 

From this result one can nicely see that, in the limit when m 7 — > 0, the electric field vanishes as 
El oc m 2 , rendering the longitudinal polarization unphysical. On the other hand, when m 7 grows 
and becomes comparable to, or larger than, k/a, the amplitude of El is unsuppressed. 
Writing st, £l as 

£{t,l) {k, rj) = c (T)L ) (k) a(r,L) {k, rj) ex P ( - * / ^(t,l) rf)dr/\ , (77) 
where a(T,L) an d ^(t,L) are rea l by construction, one obtains from ((7i|) and 1)75(1 



^ + 2^ = , ^ + k 2 + m y=u; 2 , (78) 

U>T CUT OLT 



and 



^ + 2^ ^ , (79) 

ul «l A; 2 + m 2 a 2 

cJ'r 7* o 2Hak 2 a' j 9 , , 

_L +A; 2 + m 2 a 2 + _L = ^2 ^ (g0) 

«L k 2 + m^a 2 «L 

where the prime denotes a derivative with respect to ry (' = {d/drjj). Solving the first equation of 
(|7H|) and Eq. (Z^l gives 

1 1 / p Ha'k 2 a 

«r = ^= , ai = - F =expl-/ ^- ——drj). (81) 



a'r 2iTafc 2 



y^r ^ k 2 + m 2 a' 2 

Performing the integration in the second equality gives 



1 

oil = 



m 2 + (^) 2 



(82) 
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where ao = 0(770). Thus in the oscillatory regime, where one can find real solutions for ul, 
longitudinal modes become suppressed compared to transverse modes during inflation, unless they 
are non-relativistic, which implies ||/c||/ao <C m 7 . Upon plugging (|8l|) into the second equation of 
(|78|) and Eq. (|8U|). we find the following equations for ujt = u>T(k,rj) and ujl = u>L(k,rj) 



9 - 9 9 9 3/u/r\2 Iuj'/ H 2 a 2 k 2 {m 2 a 2 -2k 2 ) 

uj 2 = k 2 + aV^ + - -M ^ + 7 . (84) 

7 4W; 2w L (k 2 + m 2 a 2 ) 2 

In adiabatic limit, when 

Eqs. (|83|) and (|84|) can be iteratively solved. The conditions (|85[) are satisfied for non-relativistic 
photons (fcph = ||fc||/a "C "i 7 ), when 

<C 7n 7 (non— relativistic limit : fc p h <C m 7 ) . (86) 

Then the last term of Eq. (|84|) can also be neglected. In the relativistic case (/c p h S> m^) we get 
from (|83|) (fcph 3> (//m 7 )2, (Hm 2 )~3), but this is not sufficient for longitudinal modes. We need 
the stronger condition 



(relativistic limit : fcph 3> m 7 ) , (87) 



^ph S> (-ffm 7 )2 (transverse modes) 1 
fcph S> -f^ (longitudinal modes) J 

for the adiabatic approximation to work, and the last term in Eq. (|84|) can be neglected. 

When the conditions (|86|) or (|87jl are met, one finds from Eqs. (jHSHSJ to leading order in 
derivatives 



W (T,L) - ^0 , = y k 2 + a 2 m 2 . (88) 

The propagation speed of massive photons is then given by the group velocity 

_ d^ph _ fcph , * 

^group — — i V ot V 

dfe ph ^ph 

which is always smaller than the speed of light in vacuo. When ||fc p h|| '-'group can be <C 1. 

Here we defined the physical frequency, u p h ~ ujq/cl, and the physical wave vector, /c pn = k/a. 

When adiabatic approximation (|85|) breaks down, one has to solve for W(t,l) = "^(t.l)^ 7 ?) 
exactly. As an example, adiabatic approximation for longitudinal relativistic modes breaks down 
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when /cph approaches 2H, even if kw^ ^> (Hm^)? , is satisfied, which means that the derivative 
terms in Eqs. (|83H84j) can be neglected. In this case, and when m 7 <C k p h,H, Eqs. (j%4"|) and (|75|) 
take the following form 

u;| ~ P - 2H 2 a 2 , (90) 
(d 2 + k 2 - 2a 2 H 2 )ae L = 0. (91) 

From (|90[) one would not expect oscillatory behavior for k p h < 2H, but solving (|91|) gives 

Ehik, v) = k~ 2 + V 2 ex p| ~~ i^kr] + arctan (— ))}) (92) 

which is oscillatory. Here we used k = \\k\\ and k = k/k. Comparing this expression to Eq. (|77|) one 
obtains ujl = k — k/(l + -rrnp) for the dispersion relation, such that the group velocity becomes 

^group — 1 7 79 ' (93) 

[i + (W#) 2 ] 

This implies that subhorizon photons (/c p h > H) propagate superluminally, while superhorizon 
photons (/c p h < H) are subluminal. When one considers propagation of relativistic photons on 
subhorizon scales, one finds that longitudinal (transverse) photons propagate superluminally (sub- 



luminally). This phenomenon is similar to birefringence 
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381 ] , where one of the tranverse polar- 



izations may propagate superluminally, and another subluminally. 

V. PHYSICAL IMPLICATIONS OF OUR RESULTS 

In this section we study the cosmological consequences of a massive photon in inflation. In the 
Introduction we argued that a large photon mass during inflation (|7U|) may dramatically influence 
the photon dynamics, perhaps the most striking being the speed of its propagation: the scalar 
vacuum fluctuations act as an 'aether', which drags photons, and consequently may dramatically 
slow down propagation of light. 



Another interesting consequence may be creation of magnetic fields on cosmological scales with 
potentially observable magnitudes 0, Q|. By following the derivation of Ref. Q|, 

1 



we arrive 

at the following estimate for the (volume-averaged) magnetic field on a (comoving) scale t c 



(t,£ c ) = f ^ a \ 4 ( <L_ S\ 2 __i — ( m t ne Gaussian system), (94) 

\7rz e q^ \27rm ( f ) c i ) t c 

Bit, Q = #(— ) * (^s) 1 HK \ + Z? (in SI) , (95) 
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where we have reinserted the physical constants. Here z eq ~ 3200 denotes the redshift at the matter- 
radiation equality, Hq ~ 2.3 x 10 -18 Hz is the present time Hubble parameter (which corresponds 
to Hq ~ 71km/s/Mpc), H ~ 10 13 GeV/h is the Hubble parameter during inflation, the fine 
structure constant a = 1/137, and /xo is the magnetic permeability of vacuum. The result (|94[) 
is derived by making the assumption that 1/2 of the energy stored in vacuum fluctuations during 
inflation is converted into the magnetic energy at the second horizon crossing during radiation era, 
when also the photon mass is assumed to vanish nonadiabatically. For an alternative derivation 
we refer to Refs. 
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4Cfl. where continuous matching of the field amplitude and its derivative at 
the inflation-radiation transition are employed, and the photon is assumed to become massless at 
the inflation-radiation transition. 

Provided the scalar field is sufficiently light, naivelyfl seems that the field strength « can 
be significantly larger than the one obtained in Ref. Jjj, where the photon dynamics coupled to 
a massless minimally coupled scalar was considered. Taking account of the more recent results 
obtained in however, which correctly treats the late time asymptotic dynamics of the photon 
field in inflation, we conclude that one would get equally strong magnetic fields from photons 
coupled to a minimally coupled massless scalar, provided inflation lasts a sufficiently long time. 



Evaluating Eq. (jMJ) with ~ 100 GeV j47|, we get 

Assuming galaxy formation at z ~ 10 and an amplification by a factor ~ 80 through field com- 
pression during the collapse of the proto-galaxy 4J]|, one gets field strengths of approximately 
B ~ 10~ G, for galactic magnetic fields after structure formation at the scales relevant for galac- 
tic dynamos. i c ~ 10 kpc has been used, which corresponds roughly to ~ 100 pc physical length 
after the field compression at z ~ 10. This is most likely sufficient to seed a dynamo mechanism, 
which generates the micro-Gauss galactic magnetic fields observed today 
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APPENDIX A: EXPANDING THE SCALAR PROPAGATOR 



We want to expand the scalar Feynman propagator 



iA(x, x') 



(4ir)-*T(§) 



2 V A H D-2 2i?i 



D-l D 1 D y 

h ZA ZA — , 1 

2 ' 2 ' 2 ' 4 



(Al) 



in the modified de Sitter length function y, s and e = 4 — D around y = 0, s = and e = 0, where 



D-l D-l 1 m z ,+£R 

V = z S = — r — r rr7; 1" O 



2 " 2 D-l H 2 
Using the following transformation formula 



2 Fi h za v, — , 1 

- i 2 1 ' 2 2 4 



# 2 



(A2) 



r(f)r(i-f) 

Y{\-u)Ya + v 



2F1 

r(f)r(f -l) 



z? - 1 D-i 



za — , — 

' 2 ' 4 



( | ) f-ir(^- J ,)r(^ + I ,) 



D-l 



,11 D y 

2 F 1 [- + v,--v,2 



2'4 " 



(A3) 



and the series expansion of the hyper geometric function 

r(c) ^r(a + n)r(6 + n)z r ' 



2F1 (a, 6, c, 2) 



r(a)r(6) 



E 

n=0 



T(c + n) 



n! 



(A4) 



one can easily derive 
iA(x, x] 



u»->- r(i-f)r(f) » r(^ + i/ + n)r(^ ->/ + ») 



(47T 



y\n 
4) 



n=0 



r(f + n) 



1 r(§ - i)r(2 - f ) g r(| + u + n)r(± - „ + n ) (f )" 



Making use of r(± + 6)r(| - 6) = r(6)r(l - 6) = ^ 

in (|A2|) one obtains 

B)) 



(A5) 



r(2-f +n) 

and inserting in the first equality 



J ff D - 2 COs(7r l 



(4vr)T sin(vr§) 



^rXD-l-s + n)r(s + n)(|) ri 



n=0 



r(f + n) 



n! 



r(f - s + n)r(i - § + s + n) (|r 



n=0 



r(2-f +n) 



7?! 



(A6) 



22 



This is still formally exact in s and e. Upon pulling the n = term out of the first sum, the n = 
and n = 1 terms out of the second sum and shifting the second sum we find, 



iA(x, x') 



H D 2 COs(-7T 



D-l 



b)) 



(4tT 



sin(7r-?) 



2 

oo 



^ 1%D - 1 - s + n)r(s + n) (f 

n=l 



r(§ + ») 



77 ! 



,y )2 _ f g£(£-s + n + l)r(2-f + s + n) (|)« 



+ 



n=l 

r(D - i - s)r(s) 



r(3 



17 



+ n) 

- £> 



(n + 1)! 



i r(f - S )r(i - f + s) 



ay 



r(2-f) 



,y,2-f £(j - s + i)r(2-f + s) 



r(3-f) 



(A7) 



Then expanding in s gives 

H 2 ~ £ rr(3-e)rl 



iA(x, x') 



(47T) 



{-+7TCOt(y) -7e-^(3- £ )+C(s,£)} 



r(2-|) 

-(J)^{i + s(.cot d ) - .0 - §) + .(f)) } + -Ljrfi - f) 



oo 

?1=1 



2 - n(n + 2) ln(f) 



+ 



77 - 1 



°(*)}(|)' 



+ 0(s 2 



(A8) 



where 7e ~ 0.577 is the Euler-Mascheroni constant, ip is defined by ip(z) = V {z)/Y{z) and C(s,e) 
is an order s term that is independent on the coordinates. The 0(e) term in the infinite sum, will 
not be needed for the regularization of the vacuum polarization. 

APPENDIX B: FEYNMAN RULES OF SCALAR QED IN POSITION SPACE 

The position space Feynman rules of scalar QED in an arbitrary D-dimensional space with 
metric g^ v are 



iA(x, x') 



' X 



i^g(x)g^(x)(d^-d^), 




-2ie 2 ^gJx)g^ v (x)5 D (x - x'), 
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v = -i5Zd p {^gJx)[g» v {x)gP°{x) - g^ (x)g^(x)]d'J D (x - x')), 
x x' 

where g = det[g^ u ], acts only on the incoming and <9° ut only on the outgoing propagator. 
These rules can be derived from the Lagrangean density (|T|). In de Sitter D-space iA(x,x') and 



g pv (x) are given by (|3fi|) and (0 , respectively, and \J — g{x) = a D . 

APPENDIX C: THE 3-POINT AND 4-POINT INTERACTION CONTRIBUTION TO 
THE ONE-LOOP VACUUM POLARIZATION 

The 3-point and 4-point (seagull) interaction contributions to the one-loop vacuum polarization 
shown in figure ^(1) and^(2) are 

i[^U u ]^(x,x') = -2ie 2 a D - 2 V ' MU iA(x,x)S D (x - x'), (CI) 

and 

i[^rT] (2) (x,x') = 2e 2 a D ~ 2 a' D - 2 t P n v<T [{d p iA{x,x')){d' (T iA{x,x')) - iA{x,x')d p d' a iA(x,x')]. (C2) 

For iA(x,x) in (|C1|) . Eq. (|4"T|) will be used and iA(x,x') is given by (jHEJ- The 0(e) terms, which 
are not needed for dimensional regularization, can be dropped because they will not contribute 
when taking D to 4 later. Dropping these and all of 0(s 2 ) terms, we find for the derivatives of the 
propagator 

d p iA(x,x') = {(3f'(y) + *-jSaW(v) + h'(y)]}d p y, (C3) 

d' a iA(x,x') ={(3f'(y) + *j^W(y) + h'(y)]}d' a y, (C4) 

d p d'jA(x,x') ={Pf"(y) + B-^tf'iy) + h"{y)]}{d p y){d' a y) 

+ {(3f'(y) + sj^[g'(y) + h\y)]}d p & a y. (C5) 

j3, /, g and h are given in Eqs. Q37H4U[) in section [nj Note that the order 1/s of f(y) is constant in 
y. Keeping only terms up to order s°, we can write the portion of (|C2|) within the bracket as 



[(d p iA(x, x'))(d'jA{x, x')) - iA(x, x')d p d'jA(x, x')] 
= [P 2 (f 2 ~ ff") ~ 2^ (9" + h")] (d p y){d>y) + [ - (3 2 ff - 2j^(g' + h')}d p d' a y. (C6) 
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Differentiating the modified de Sitter length function y gives, 

d p y =aH[y5° p + 2a'HAx p + 2ia'H sign (77 - j/)<5J<5], (C7) 

d' a y =a'H[y5° a - 2aHAx a - 2iaH sign(r ? - rf)$8\ , (C8) 
d p d' a y =aa'H 2 [y5 p 5° a - 2aHAx a 5° p + 2a'HAx p 5° - 2 Vp(T 

- 2iaa'H 2 \r l - r/\S%%6 - 4i<% - ri')5° p 5°J} , (C9) 

where Ax p = x p — x' p . When taking 5 to zero only the final order 5 term contributes, and only 
when it multiplies /' 

-4iaa'H 2 8{ri - r,') lim f'(y)8 = -^^^{x - x'). (CIO) 

This contribution subtracts off the temporal component of I)C1|) . Dropping O(s) terms and insert- 
ing (|41|1 into (|(J1() . we find for the remaining part i^IF]^) 

iCn^x, a:') = -2ie 2 a D - 2 pt U K5 D (x - x'), (Cll) 



where K is the constant part of f(y), 



R _ 2^r(3- £ ) 



1 / 7T8 \ . . 

- + 7rcot (—J - 7e - ip{3 - e) 



i + J--| + ln(2) + 0( e ). (C12) 
e zs 4 



4T2(2-| 

We shall call the left over portion of (|C2|) . i^n^]^, x') and use a bar over the metric tensor 
to indicate that its zeroth components have been removed, ff v = rf" v + 5q 5q . Calculating the 
derivatives of f(y) and dropping some unimportant 0(e) terms one easily finds 

e y 4-e £ yi-e y 

and 



f/ 2_ ff „__J: 1 (2-§)JT Wl lnr ^ 1 1 



//= "1^|^ + 2 £ (1-|) ^-^f + ^4 ln (4)-4l + 2J- (CM) 
Differentiating Q39J1 and keeping only terms up to order e° yields 

g'(y) = - , g"(y) = -^- (ci5) 
y y 2 

When we perform the summation in (|4Uj) and calculate the derivatives of the resulting expression, 
we obtain 

y-4 + 2(y-6)ln(g) 

^ (y) = ^734)2 — -> (°i 6 ) 

, 48-16j / + 16j/ln(|)+^-2j / 2 ln(|) 

* {y) = y¥^F • ^ 
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Upon inserting into and (|C13HC17|) into l|C6|). one gets 



i["Il£](x,x') =2e 2 a 3 - e a' 3 - e ^ 2 /3 2 



Aaa'H 2 Ax^Ax u 



I_|y4- B e y 3 ~ e y 3 -i y 2 U M ; 
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1 1 (l- £ )(2-|) 1 



1 - f y 3 ~ e 



+ 



K 



1 /l 



1 3 
4s~ + 4 

'2/ 



1 1 



1 
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where here and in the rest of the appendices upper indices are raised by the Minkowski metric r]^ u . 
The portion of (|C18j) within the first bracket and the constant in front of it can be written as 
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We want to write the vacuum polarization in a form which makes explicit that d fJi i[^H u ](x, x') = 
d' I/ i[ fl Il' / ](x , x') = 0. To find such a manifestly transverse form we will use the following identities 
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When inserting ()C22|) into (|C19|) . the local term exactly cancels i[ )1 Il'(](x, x'). We find for the 
remaining part of (|(J19|) . 
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where a = e 2 /47r is the fine structure constant, and we dropped all 0(e) terms, except in the first 
line. These terms are not required for regularization. Note that the second term is not yet in 
manifestly transverse form, because of the factor standing to the left of the projector operator. 



Bringing it to the right and combining with the rest of (|C18|) we get for the sum of all terms in 
Eqs. (EU) and (IH2l) 
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(y-4)3 

We will use the following ansatz for the terms in (|C24|) . which are not in manifestly transverse 
form 
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Comparing this to the relevant terms in (IC24I) . one finds 
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One can find the following solutions to these differential equations 
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where the dilogarithm function, defined by Lii^z) = — Jq ^j-^dt appears. Upon inserting these 
solutions into ()C25|) and combining the result with the first two terms of (|C24|) . one obtains the 
following expression for the graphs in Fig^(l) and Fig^(2) 
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APPENDIX D: INTEGRATION OF THE RETARDED VACUUM POLARIZATION 

AGAINST THE PHOTON FIELD 



In order to simplify the effective field equation for the photon field, we need to evaluate the 
integral 
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We can simplify the integral by using the identity, 
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where / can be an arbitrary function of x — x' and a bar over a tensor is used to indicate that its 
zero components have been removed, e.g. fj^ = n^ u + 5q<5q, B^ = 8^ — 5qB°. Using this we find 
for the integral (|D1|) 
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In the last term the transverse projector [ ^P^] hits a function of x — x 1 only. Therefore it can be 
replaced by -(ry^V 2 - B»B V ). Then this expression can be recast as 
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One can get rid of the step functions by choosing appropriate limits of the integration. Then one 
can rewrite (ID8I) as 
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where 6 is the angle between Ax and k, r = \\Ax\\ and k = \\k\\. Performing the r and cos# 
integration yields 
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From this expression we obtain 
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Using (jDllj) one can write the contributions to the \i = component of Eq. ()D7|) as 
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Upon using (|D4|) and performing partial integration this can be recast as 
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When performing the (% differentiation all integrals cancel, only one term remains. It comes from 
the differentiation with respect to the upper limit of the integral drj'eo (k , rf) ^ cos (^ At? ) ; thus 
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The terms contributing to the spatial components of (|D7|) can be written as 
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where ([DIOjl and (jDllj) have been used. We now insert Eq. (|D4|) into the terms in which ik u e u (k, rf) 
appears and partially integrate. When we then perform the remaining differentiations all integrals 
cancel, and only a term, which comes from the differentiation with respect to the upper limit of 
the first integral, remains. This term equals 
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Combining this with ()D13|) we finally find for the integral ()D1|) 
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